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1 abstract 


Some of the potential advantages of fiber reinforced composite materials include: high 
strength at elevated temperatures, high strength to weight ratio, improved toughness, 
improved creep and fatigue strength and controlled expansion and conductivity. They 
also offer the unique advantage to tailor the properties and shape of the high perfor- 
mance material for the application desired. Therefore, there is a critical need of fibers 
capable of effectively reinforcing these intermetallic matrices at elevated temperature. 
This research work describes the CVD synthesis of such libers which can be used for the 
potential applications in high temperature composite materials. This process was used 
due to its advantage over other fiber synthesis processes. It is extremely important to 
produce these fibers with good reproducible and controlled growth rates. However, the 
complex interplay of mass and energy transfer, blended with the fluid dynamics makes 
this a formidable task. 

The design and development of CVD reactor assembly and system to synthesize 
TiB 2 , CrBi , B 4 C , TiC fibers was carried out during the course of this research. The 
thermodynamic and kinetic analysis of T 1 B 2 and CrB 2 , was done to locate the windows 
for the pure deposit and to scan out different precursors by keeping the emphasis on 
effect of temperature, concentration of reactants, gas flow rate and morphology. Residual 
stress analysis for estimating stresses arising from thermal expansion mismatch has been 
determined. Various techniques to improve the mechanical properties were also carried 
out. The developed fibers showed extremely good strength. The TiC, B+C, and B fibers 
matched the properties of best known results. TiB 2 with the intermediate spongy layer 
showed the strength up to 250KSI. The exploratory work on CrB 2 gave excellent result 
with initial strength values varied from 60 to 170 KSI. Various techniques for improving 
the fiber properties has been elaborated. The crystal structure and its orientation for 
TiB 2 fiber is discussed in. In short, this report gives an overall view of CVD process to 
develop CrB 2 , TiB 2) and other high performance ceramic fibers. 
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2 Technical Background 


A new generation of composite materials is revolutionizing today’s aircraft and auto- 
motive industries [Edie, 1987]. In these high technology applications, weight, strength, 
and stiffness are critical parameters which can be achieved by using high-performance 
ceramic material composites [Edie and Dunham, 1987]. These ceramic materials offer 
many outstanding properties for high- temperature applications. This attention led to 
the discovery of fibers which have excellent high-temperature strength and corrosion 
resistance, good thermal shock resistance, low density, and low coefficient of thermal 
expansion. High-performance fibers are available only at high prices and, as a result, 
tend to serve specialized and restricted markets. Basically, three factors affect this price: 
raw material, manufacturing, and more important new market development cost. Sec- 
ondly, the use is partly limited because of their brittle nature and large degree of scatter 
in strengths which often lead to catastrophic failure. A number of techniques have been 
considered to lessen or eliminate these problems [Richardson, 1984]. 

Over the past few years, numerous patents and publications have revealed the 
use of simple organo-metallic compounds, metal halides, and polymer precursors for 
routes to ceramics fibers [Campbell et al., 1949]. High Gibb’s free energy of formation 
of organometallic compound frequently allows for lower temperature of formation of 
ceramics. This fact is utilized in the CVD process. The volatile organometallic precur- 
sors allow for simple purification by superfractionalization so that CVD process results 
in an ultraclean material. Very recent experimental observations predict that traces 
of impurities in a fiber (e.g., calcium) result in a fast high-temperature degradation 
(Johnson et al., 1987). 

The achievement of consistently high fiber strengths over a long length is impor- 
tant for successful implementation in many composite systems. This requires highly 
controlled processing operation and conditions [Brenman and Prewo, 1982], Coating 
defects may also induce some kind of flaws or microbending losses in the fiber [Glogo, 
1975; Gardner, 1975]. Only two ceramic fibers are available commercially (boron, silicon 
carbide). But these fibers are not matching the properties required for the advanced 
composite system. For the intermetallic matrix reinforcement we are looking for a spe- 
cial kind of fiber. Because metal matrix composites offer low weight, high strength, high 
stiffness, and high temperature capabilities to satisfy the performance requirements of 
future advanced aerospace -systems. By choosing a matrix material that bonds well 
to the fiber, is ductile, and stands up to the environment, properties of the compos- 
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ite can be tailored to the needs of advanced power and propulsion systems. Titanium 
diboride, titanium carbide, boron carbide, and chromium diboride fiber compositions 
were synthesized as reinforcement material for intermetallic matrices in this research. 


3 Objective of the Work 

Fiber reinforced composites occupy a central role in the development of new materials 
and has overcome many of the limitations of traditional materials. Some of the potential 
advantages of fiber reinforced composite materials include: high strength at elevated 
temperatures, high strength to weight ratio, improved toughness, improved creep and 
fatigue strength and controlled expansion and conductivity. They also offer the unique 
advantage to tailor the properties and shape of the high performance material for the 
application desired. 

Therefore, there is a critical need of fibers capable of effectively reinforcing these 
intermetallic matrices at elevated temperature. This report describes the CVD synthe- 
sis of these fibers which can be used for the potential applications in high temperature 
composite materials. The CVD method has several advantages over other conventional 
methods; hence, it has been adopted for this research. It is extremely important to pro- 
duce these fibers with good reproducible and controlled growth rates. Other properties 
axe to be controlled include thickness, composition, purity, crystallinity and surface 
morphology. These properties are highly dependent on the deposition condition. Hence 
it is necessary to correlate the fiber properties with temperature of the substrate, sys- 
tem pressure, gas flow rate and concentration of the reactants. However, the complex 
interplay of mass and energy transfer, blended with the fluid dynamics makes this a 
formidable task. 

Hence, goals were to design and develop the CVD reactor systems which can syn- 
thesize TiB 2 , Cr B 2 , B^C, TiC fibers. The thermodynamic and kinetic analysis of TiB 2 
and CtB 2 , with the emphasis on effect of temperature, concentration of reactants, gas 
flow rate and morphology were carried out. Investigation on the mechanical properties 
of the developed fiber including strength and modulus were conducted. Residual stress 
analysis for estimating stresses axising from thermal expansion mismatch has been de- 
termined. Various techniques to improve the mechanical properties were also carried 
out. 
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4 CVD Synthesis 


Broadly defined, chemical vapor deposition, or CVD, is the formation of solid products 
via chemical reactions of gaseous precursors. A typical CVD process would involve a 
dynamic flow system in which gaseous reactants pass over a heated substrate. The 
gases react chemically to produce a condensed coating on the substrate plus product 
gases, which, together with any remaining reactant gases, dynamically exit from the 
hot reaction zone. 

The application of CVD processes centers around the production of thin films for 
semi-conductors and other solid state electronic devices (Hess et al., 1985), coating 
of cutting tools and surfaces needing erosion and/or corrosion protection (Yee, 1978), 
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coating of fibers used in forming composite materials (DiCarlo, 1985 and Guinn and 
Middleman, 1989), and containment coatings for nuclear fuel and nuclear waste particles 
(Spear, 1982). CVD processes are also used in the manufacturing of objects with com- 
plex shapes (e.g., refractory crucible) out of materials such as tungsten, molybdenum 
and rhenium which resist conventional machining and fabrication techniques (Iwasa et 
ah, 1987; Shinko and Lennartz, 1987). 

A vast number of chemical compounds are considered for these applications (Blocher 
et ah, 1984). Important materials being prepared by CVD for electronic and optic de- 
vices include Si, Ge, SiO 2 , GaAs, CdS, ZnSe and related compounds. Supercon- 

ducting materials such as Nb 3 Sn 7 NbzGe and NbC x N y can be deposited as well (Wahl 
and Schmanderer, 1989). Erosion and corrosion resistant coatings commonly applied by 
CVD to tooling materials are TiC, TiN and Al 2 0z- In addition, an appreciable amount 
of CVD research is being performed on coatings of TaC, TaN, W 2 C , TiB 2 , SiC, Si z N A:i 
B 4 C, BN and B for use on cutting tools, tubes, fibers and applications in which the 
materials are subjected to high temperatures, and corrosive and erosive environments. 


4.1 Synthesis of High Performance Ceramic Fibers by CVD 

The CVD process involves the formation of solid by decomposition or reduction of one 
or more gaseous components upon a heated substrate. Typically, a small- diameter sub- 
strate wire is run through the glass reaction tube, and suitable gases are introduced. 
The substrate is resistance heated causing the gas to react and deposit on the heating 
wire. Uniformity of the coating mainly depends on the uniform surface structure of 
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the base fiber, temperature and rate of supply of fresh reacting gas to the fiber sur- 
face [Pierson et al., 1978]. The non-uniform temperature distribution results in heavier 
deposits on the hotter region and thereby introduce some flaws in the material. Deposi- 
tion non- uniformities due to the fluctuations in the concentration of a coating gaseous 
reaction mixture along the fiber might be overcome by lowering pressure in the system 
[Sherman, 1987]. 

The coating adhesion is excellent when some kind of physical/chemical bond is 
formed by diffusion, or reaction between the coating material and the base. However, 
this diffusion (or interaction) should be at the ^optimum” level to maintain the strength 
and modulus of the fiber. Cleanliness of the surface of the base material is very im- 
portant where the surface and base material interaction is necessary to obtain good 

t 

adhesion. Cleanliness of the matrix is not so important when only physical adhesion of 
the material is required. 

The supersaturation and temperature of the reaction process are important factors 
affecting the type of deposit. Low temperature or higher supersaturation of gases leads 
to a gas phase nucleation, thereby decreasing the efficiency of coating and adhesion. 
Higher temperature or low supersaturation gives epitaxial growth or dendrites kind of 
growth which is determental for the coating process. Hence, optimum temperature and 
supersaturation is necessary to obtain the required deposit such as amorphous deposit 
or fine grain polycrystal [Naslain et al., 1979; Schlichting, 1980], 

Deposition rates vary with the material deposited and the process used. Coating 
efficiency is the measure of material deposited to material supplied as volatilized com- 
pound, might vary greatly with the coating process and strongly depends on operating 
conditions and geometry of a reactor. 

Generally, the material necessary for the CVD process can be supplied by a source 
gas. For the CVD process, two general reaction routes can be followed: 

• Reduction of volatile halides 

• Pyrolytic methods 

The hydrogen reduction of halides process is an effective technique for the synthesis of 
ceramic fibers by CVD. The principle of this technique is as follows: Hydrogen saturated 
with metal halides (e.g., TiCU , ZrCl 4 , etc.) or nonmetal halides (e.g., CCl A , BCl 3 , 
etc.) or mixture of these halides is passed into a reactor in which the base filament is 
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heated (resistively or by any other method). Typically, for example, for a TiB 2 filament 
synthesis, the following reaction occurs: 

TiCl 2 + 2 BCh + 5 H 2 ^TiB 2 + 10HCI (1) 

The pyrolytic method represents a thermal decomposition of an appropriate com- 
ponent with low boiling point. For example, silicon carbide fiber can by synthesized 
from CHz.SiClz 


CH z SiCl z -*SiC + ZHCI (2) 

The surface deposition reaction rate varies with the material deposited, the reactant 
composition and temperature. Usually a material is deposited with the rate of 5 to 100 
/zm/hr. The important parameters which should be considered in this process are: 

• Purity of the deposit and uniformity of coating. 

• Coating adhesion. 

• The effect of supersaturation and temperature on the structure. 

• Efficiency of CVD reactions. 

• The reactor geometry and operational arrangement. 


5 Experimental Set-up and Operations: 

The various applications of CVD have resulted in the development of various kinds of 
reactors, heating systems and feeding systems. The manner in which the substrate is 
heated is the primary factor in the design of the deposition chamber. The design will 
be influenced by the geometry, shape, and the composition of the substrate, the type 
of deposition process used, the nature of the coating desired and by the economic and 
persona] preference factors. The most important criteria among all these should be 
the performance and cost. The CVD system used is illustrated in Figure 1. The main 
components of the system are: the CVD chamber, the vapor precursor feed system, the 
temperature monitoring unit and the effluent gas handling system. However, the CVD 


12 




Figure 1: CVD System Lay-out jj 

reactor is the most important part of the process. Our batch system basically consists 
of a quartz reactor of 1.5 inch diameter with projections for electrode insertion. The 
distance between the electrode is 15 inch (i.e. sample fiber length). The electrodes are 
of cylindrical shape with a small hole to mount the substrate fiber. The bulkiness of 
the electrodes was necessary to avoid overheating of the electrodes. The extra heating 
coil was wound around the reactor to avoid the condensation of sub- metal chlorides on 
the walls of the reactor. The typical wall temperature was varied from 200-600°C. The 
reactor is shown in Figure 2. 

The vaporizer was maintained at 250°C. This also serves as a mixing chamber 
for the reactants and the carrier gases, which is critical for uniform deposit. The 
temperature of the substrate was measured by an optical thermometer (±5°C accuracy) 
and controlled by monitoring the applied voltage. The de-oxo system contains activated 
copper turnings maintained at 500°C to remove traces of oxygen from the carrier gases. 
The traces of moisture were eliminated by passing through Dehydration unit (anhydrous 
calcium oxide). The liquid reactants were pumped through a graduated peristatic pump. 

In a typical operation, j;he fiber was mounted between the two electrodes and sys- 
tem was made air-tight. The reactor surface and preheater were brought to the required 
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temperature. Meanwhile, the reactor was continuously flushed with extra pure argon 
gas which was pre-purified by passing through de-oxo and dehydrating units. After 
making sure that there is no more oxygen inside, the hydrogen flow was started at 1.5 
lit/min., again passing through the pre-purifying unit. The purity of reactants and 
carrier gases is extremely important because the trace amount of impurities has a dele- 
terious effect on the physico-chemical properties of the fiber [Powell et. al., 1955]. The 
fiber was brought to required temperature (1000-1300°C) under hydrogen atmosphere. 
Once the temperature was stabilized, reactants were pumped into the vaporizer/mixing 
chamber, for a fixed amount of time. Once the deposition was over, the fiber had been 
cooled to room temperature slowly under hydrogen atmosphere to avoid the thermal 
stresses. 


5.1 CVD - continuous Process 

The synthesis of ceramic fibers by CVD technique will be carried out in continuous 
manner. The set-up is shown in Figure 3. 

The metallic or graphite filament to be coated enters the apparatus through a 
mercury seal and is heated in an inert gas atmosphere in the inlet chamber. The 
hot fiber enters the main plating chamber through a small hole in the inlet chamber 
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and finally leaves the apparatus through a similar arrangement at the exit end. The 
procedure is the same as that of batch operation. 

The wire is heated resistively. The reacting gaseous mixture enters the reactor at 
different axial and angular positions to increase accessibility and uniformity of reactants 
on the filament surface. A motor- driven reel draws the fiber through the apparatus at 
a steady, controlled rate. 
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6 Materials and Chemicals 


6.1 Substrate selection: 

The final properties of the CVD fiber are dependent on the substrate fiber used. Several 
substrate uses are reported in the literature [Person and Randich, 1978; Pierson et. al., 
1979; Bryant and Meier, 1974]. Basic criteria for a satisfactory substrate selection are 
as follows: 

1. Ability to withstand the temperature of deposition. 

i 

2. Chemical inertness to the deposition atmosphere. 

3. Thermal expansion is close to that of the deposit. 

Several materials were tested for our substrate. Among them Ti, W, Mo, Cu, Ni 
and C. showed good potential. In all of the experiments, hydrogen is the carrier gas. 
At the deposition temperature, all the material except W and C, absorbs hydrogen to 
various degree to form a metal hydride leading to embrittlement. Even though most 
of the above mentioned material are resistant to HC1, only W and C can withstand 
prolong attack at deposition temperature. Finally, ductile materials like Ni, Cu are 
having high thermal expansion mismatches with the ceramic deposit compared to W, 
Mo and C, creating a larger residual stresses in the deposit. After analyzing all these 
results, we selected W, C, Mo as the substrate materials for the synthesis of fibers. The 
properties of the commercial fibers/ wires used are given below. 

C (Thornel, From AMCO) Fibers: 

• - PAN based: Dia.= 5-10 /xm, Surface area = 0.5 m 2 /g, Strength=100-300 KSI, 
2000-5000 filaments in a bundle. 

• - Pitch based: Dia.= 6-10 ^m, Surface area = 0.3 m 2 /g, Strength =300-400 KSI, 
2000-5000 filaments in a bundle. 

W (wire Johnson Matthey and GE): 

• Johnson-Matthey fiber - dia. = 100 /xm, tolerance ±8%, Surface area = 0.07 m 2 / g 
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Table 1: Reactant Analysis 


Reactants 

Impurities in parts per million(PPM) 

H 2 

0 2 = 0.5, Cl 2 = 0.1 

Ar 

H 2 — 0.1, O 2 =0.1 

BCl 3 

Si=8, C7 2 =80, COCl 2 =200, Al=30, Fe=22 

TiCU 

Fe <10, Si=50, Sn=100, C=50, V=10 

ecu 

S=10, H 2 0 = 20 


• Johnson-Matthey fiber - dia. = 75 fim, tolerance ±8%, Surface area = 0.07 m? j g. 

• GE fiber, Electro polished and Electro etched (customized), 

“ Electro etched - Dia. = 40 /xm; tolerance zL3%, Surface area = 0.5 m 2 /g 

— - Electro etched - Dia. = 100 /x m, tolerance ±3%, Surface area = 0.5 m 2 /g 

— - Electro polished -Dia.= 70 /x m, tolerance i 3%, Surface area =0.3 m 2 /g 
“ - Electro polished - Dia. = 70 /xm and 100 /xm, tolerance dh 8% 

Ti (Wire, Johnson Matthey): Dia. = 100 /xm, tolerance ± 8% 


6.2 Chemicals: 

The quality of the film is directly dependent on the purity of the reagents used. It is 
of utmost importance in CVD to use extra pure material [Powell et ad., 1955] to avoid 
poor physico-chemical properties to the film. Following chemicals are used: 

• Hi and Ar obtained from CS Gas Products, Linde Div. of Union Carbide (elec- 
tronic grade, extra dry, 99.99% pure). The gases were further purified by passing 
through de-oxo and dehydration units. 

• TiClt , BCl^y CrOiCliZLiid. CCl A obtained from Johnson Matthey Co. are of 

commercial grades. Wherever possible, these materials were redistilled to improve 
the purity. The typical chemical analysis is given in Table 1. 
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6.3 Pretreatment of the Substrate 


One of the main requirements of the deposited films is that it should adhere to the core 
at all operating temperatures and capable of withstanding severe deformation without 
spalling. Hence surface cleanliness of the substrate is of utmost importance. This was 
done three different ways: 

1. treatment with boiling concentrated nitric acid for 10 mins. 

2. treatment with steam at 250° for 10 mins. 

3. treatment with H 2 at 500°C for 5 mins. 

These treatments riot only cleans the surface of the substrate material but also 
increases the surface area. The best results were obtained with hydrogen treatment 
where surface area of W was increased by 10% and that of carbon was by 25%. This 
method was used in most of our sample runs; if otherwise, mentioned 


T Titanium Diboride Fiber Synthesis 

There is a critical need of fibers capable of effectively reinforcing ceramic and metal 
matrices at high temperature. The high specific strength and stiffness at elevated tem- 
peratures combined with improved toughness, creep and fatigue resistance of these com- 
posites, overcome some of the limitations of traditional materials. Titanium diboride 
(TiB 2 ) is the most promising candidate for reinforcing the metal matrix. It processes 
high strength at elevated temperatures, and exhibits high stability and hardness. Its 
compatibility with Ni-Al and Fe-Al matrix was reported in the recent NASA report 
[Misra, 1988a, 1988b]. This research gives an overall view of an extensive work carried 
out in the Laboratory for Ceramic and Reaction Engineering, on the application of 
Chemical Vapor Deposition technique in the development of TiB 2 fiber. Detailed work 
on growth and characterization of fiber including simple kinetic expression is discussed. 
The residual stresses developed due to deposition condition are also explained. 

Even though several techniques are available for manufacturing TiB 2 fiber, Chem- 
ical Vapor Deposition (CVD) is used because of its flexibility and ability to deposit 
dense, pure refractory coating. Again, one can obtain the fiber well below the melting 
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point of the material, with good control over the morphology. Some of the commonly 
used precursors for the deposition of TiB 2 include Ti + B [Rauch et ah, 1968; DiCarlo, 
1985; Bunsell, 1988]; TiCU + BCU [Lennon, 1965; Anon, 1966, Chernyshova et al., 
1985; Hess et al., 1985; Biswas, 1986]; TiCU + BBr 2 [Graham, 1971; Chamberlin and 
Skarman, 1966]; TiCU + B 2 H& [Viguie and Spitz; 1975] and Ti^BH^)^ [Monney and 
Radding, 1982]. In short, one can summarize the precursors to three different categories 
to produce TiB 2 fibers. 

• gas-solid reaction 

• bimolecular decomposition and 

• decomposition of a organometallic compound. 

In gas-solid reaction, boron is directly deposited on titanium metal disc or a wire 
[Warth, 1923; Thebault et al., 1976] by reducing BCU . They have shown the Ti-B 
reaction to be a diffusion controlled governed by the equation x = Ct* where x is 
the thickness of the boride film and t is the time with C as the rate constant. The 
deposition of Ti from reduction of TiCU on boron was reported by Bouix et al., (1986) 
and Vincent et al, (1986). This reaction also showed that the thickness of the boride 
layer is proportional to £2 . The interaction of Ti metal with B powder is discussed 
by Neronov et al., (1981). The main disadvantage of this process is the formation of 
multiple borides along with TiB 2 . The property of final fiber is very poor with many 
defects. 

Organometallic compounds generally are borohydrides [Campbell, 1956] which are 
unstable, pyrophoric and can be easily decomposed at a relatively low temperature. The 
best precursor is Ti(BH ^) 3 which is prepared by reacting TiCU and LiBH$ [Powell et 
al., 1955]. Another precursor is titanium monochloro borohydride which is prepared 
from TiCU and Al{BH±)z which has the same volatility as Al{BH±)z. The main dis- 
advantages of this process are: 

• purity of the precursor, 

• co-deposition of free boron and non stoichiometric deposit , 

• compound instability, 

• less control over the morphology of the deposit. 
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Finally, TiB 2 is deposited on a hot substrate through a bimolecular decomposition 
of volatile compounds of titanium and boron. This process offers faster deposition rate 
and greater control on the composition of the deposit. Moers (1931) studied extensively 
the decomposition of halides of boron and titanium on hot substrate. The titanium 
precursor is always TiC/ 4 (b.p 136. 4°C) whereas boron precursor is either boron halide 
BBr 3 or Biz) or boron hydride (B 2 H G) etc.). Deposition from diborane requires 
low temperature (600-900°C) [Randich and Gerlach, 1981] but chlorine contamination 
of the deposit is difficult to avoid. The co- deposition of free boron is also competing. 
The biggest negative factor is its high toxicity and extreme flammability. Even though 
BBr$ or BI Z decompose at low temperature, it is advantageous to use metal and boron 
compounds containing a common atom. It is reported that the film deposit from BBr z 
is of poor quality [Zeman et aL, 1982], 

The most extensively studied reaction for TiB A deposition is between TiCl A and 
BCl$ Peshev and Niemyski, 1965]. Their study showed the different morphological 
forms of TiB 2 , depending on the deposition temperature. Gannon et al., [1963] ob- 
tained the pure dense deposit of TiB 2 at 1600°C and elaborated the crystal orientation 
mechanism. The use of different analytic technique was reported by Gebhardt and 
Gee [1965] and Pierson and Mullendore [1981], along with the optimum deposition 
conditions. Bessman and Spear [1975, 1977] in a series of three papers published the 
thermodynamic analysis of TiCl 4 — BClz system and also the effect of gas phase com- 
position on the morphology of TiB 2 . From our analysis, it was concluded that the best 
precursors for TiB 2 deposition are BClz and TiCl 4 . These are available commercially in 
pure form, and can be handled without any difficulty. The most important factor is the 
operatability of the system with wide ranging conditions. Most of the experiments were 
carried out on tungsten substrate because of its resistance to HC1 attack and inertness 
to hydrogen at deposition temperature. Electro polished and electro etched tungsten 
fibers were supplied by General Electric Co. with tolerance level ±3%. The commercial 
tungsten filaments of Johnson Mattheys are also tried. 

7.1 Results and Discussion on TiB 2 Fibers 

7.1.1 Thermodynamic Study 

On actual practice, prior to deposition with the various precursors, it is beneficial to 
conduct a thermodynamic -analysis to get an insight of various conditions available for 
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an experimental operation. This will reduce the amount of experimental data required 
to optimize the process. As discussed earlier, Bessman and Spear [1977] have analyzed 
the thermodynamic data of B + Ti system to obtain the window to get the pure 
TiB 2 . Randich and Gerlach [1981] found that the pure TiB 2 phase is possible only 
at high B to Ti ratio. Our operation is similar to above analysis but covered whole 
spectrum of parameters. This operation gives the composition of a reaction mixture 
as a function of operating parameters such as temperature, pressure, and composition 
assuming equilibrium established in the operating system. Program, NASAF obtained 
from NASA was used in determining the equilibrium composition of the gas mixture 
[Gordon and McBride, 1971]. The constants are calculated from thermodynamic data 
available from JANAF tables [1983]. 

Temperature, pressure and the concentration of the reactant mixture [TiCU + 
BCl z + Ar + H 2 ] were input to the calculations. The criteria of conservation of each 
element were imposed in calculating the equilibrium composition that had the minimum 
Gibbs free energy. The various species considered in the study are B, BC1 y BCl 2) BCl z , 
BH, BH 2) BCH z , jB 2 , Cl ,C7 2 , H, H 2) HC1, TiCl 2) TiCl Zi TiCU , B(s), Ti(s), TiB(s), 
TiB 2 (s), Ti 2 B$($) and Ti z B A (s). 

Figure 4 shows the phase boundary for pure TiB 2 . Here BCl z /TiCU ratio is 
plotted against the operating temperature for various H 2 jTiCl A ratios. Below the 
plotted line, one can exclusively get pure TiB 2 where as above it is the co-deposition 
of TiB 2 + B. This region of pure TiB 2 increases with a decrease in H 2 to TiCl$ ratio. 
Under our operating conditions, only possible solid phases are TiB 2 and B. Multiple 
borides (Ti 3 S 4 , Ti 2 B$ and TiB) are observed only at high temperature which is not 
the scope of this study. This plot gives the thermodynamically required concentration 
of reactants to get pure condensed TiB 2 phase. 

Analysis was carried to find out the effect of H 2 /TiCl A and BCl z /TiCU ratio 
on the efficiency of deposition at various temperature. Deposition efficiency is defined 
as the deposition rate of TiB 2 at equilibrium divided by the deposition rate in the 
absence of any thermodynamic and kinetic limitation. The deposition of TiB 2 efficiency 
increases with temperature and decreases with BCl z /TiCU ratio. However, for high 
values of H 2 /TiCl A ratio (> 20 ), the deposition efficiency seems to stabilize. Our analysis 
concludes that the optimum H 2 to TiCl± ratio is around 25 and BCl z /TiCl 4 ratio is 2 
at operating temperature [1400-1700 K]. 

The effect of pressure- on the deposition efficiency at various temperature for the 
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optimum ratio of reactants is given in Figure 5. It shows that the deposition efficiency 
increases with decrease in pressure and increase in temperature. This effect is pro- 
nounced in the lower temperature range whereas at our operating temperature range 
[1400-1700 K], it is having negligible effect. We maintained 1 atmosphere because of 
the difficulty in handling low pressure operations and the marginal effect of pressure on 
deposition efficiency. 

The typical concentration profile of reactants and product gases is given in figure 
6. at various temperatures. 

At any given temperature H 2 , HC1, BCh and TiCU make up 95% of the gasphase. 
Other gases, BC1, BCI 2 , BH 2: BH $ , TiCl 2 , TiCl 2 . H and Cl are in negligible amount. 

This thermodynamic study gave the composition of various compounds and de- 
position rate at equilibrium, however, their actual use in the experimental operation is 
limited. In general, CVD systems are often operated far from equilibrium conditions. 
This leads to a different experimental result than the thermodynamical^ predicted. So 
we decided to find out experimentally the kinetics of deposition. 
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7.1.2 Kinetic Study 


The chemical vapor deposition process is mainly controlled by kinetic, diffusion and 
thermodynamic factors. Under real conditions, these factors manifest themselves si- 
multaneously and determine the feasibility and direction of the deposition process as a 
whole. The various parameters studied include total gas flow rate, temperature of the 
substrate and concentration of the reactants. The operating pressure was maintained 
at 1 atm. 

The rate of growth can be measured by the change of weight, provided the surface 
area available for deposition is fixed. Therefore, deposition rate must be defined as 
the weight or volume change per unit surface area. Hence overall deposition rate is 
expressed as r 2 - r x /t where r x and r 2 are initial and final radius of the fiber and t is 
the deposition time. The plot of radius versus deposition time at 1250°C showed the 
linear dependence of growth time (Figure 7). 

The effect of total gas flow rate on deposition rate was studied to eliminate the mass 
transfer limitations. The plot is shown in Figure 8, where the operating temperature 
was 1250°C with H 2 /TiCl± and BChfTiGU ratios are 25 and 2 respectively. It is 
observed that the deposition rate is independent of total gas flow rate. That means 
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the rate data is solely influenced by the kinetic and thermodynamic factors and mass 
transfer limitations are eliminated. This result is similar to that of Mamet’ev et al., 
(1974). However, at higher temperature and lower total gas flow, the deposition rate 
drop indicating the change in deposition mechanism. 


The rate of deposition of TiB 2 can be expressed as a product of concentration of 
reactants and temperature, which are independent of each other. 


k = k’f(c) 

where k is the deposition rate, 
k’ is a function of temperature, and 

/(c) = Cbci 3 Cricu r where m and n are the order of the reaction with respect to BClz 
and TiCl A respectively. 


Since the hydrogen concentration is always in excess, the deposition is assumed 
to be independent of it. Similarly the concentration of HC1, a reaction product. Since 
deposition rate is small and hydrogen is used in large excess, the net gas flow rate more 
or less remains constant. So we varied the input ratios of H 2 /TiCl 4 and BCI 3 /T 1 CI 4 
to find out its effect on deposition rate. 

The effect of BCl 3 concentration on the deposition rate is shown in figure 9 at 


25 




ln(C BCl3 ) 


Figure 9 : Effect of BCI3 Concentration on Deposition Rate at 1250 °C and 
H 2 /TiCl 4 = 22 

■■■»» ■■ ■ ■ 2 z Z "v ■ ■ ■ z • " — ; gsassa a -- as agnw pi " 1 ™ ,l 

fixed H2/TiCU ratio (22) and temperature ( 1250 °C). 

The graph is linear and slope gives the reaction order with respect to BCl$ (i.e., 
m =0.8). Similarly, the concentration of TiCl 4 is varied, keeping the BClz/TiCU ratio 
constant. This is shown in figure 10 . From the plot, it can seen that the concentrations 
of both BClz and TiCl A varied simultaneously and the slope gives the overall order (1.0). 
This gives the order with respect to TiCl A to be 0.2 (= n). We have also observed the 
decrease in deposition efficiency with increase in H2/{TiCU + BCl$) ratio which is 
similar to that of Peshev at al., ( 1965 ). 

Temperature plays an important role on the rate of deposition a s well as on the 
morphology of the deposit. Increasing the temperature, increases the deposition rate 
rapidly. This is delineated in the Arrhenius plot (figure 11) where ij (deposition rate) 
is plotted against inverse temperature readings between 1200 - 1450 °C for two different 
concentrations of reactants. It is observed that decreasing the concentration of BClz 
and TiCl 4 decreases the deposition rate but the slope ^ remain unaffected. 

Here E is the activation energy (K cal/mol k) and R is the gas constant. This 
confirms that our data is .riot governed by mass transfer limitations. The activation 
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energy obtained from above plot is 40.7 k cal/mol K. This is in an excellent agreement 
with the values obtained by Bessmann et al. [1977] [40 ±12 K cal/mol K] with the 
reduction of BCI3 and TiCl 4 by hydrogen. There is no change in slope between 1200- 
1450°C; shows that only the reaction kinetics controls the deposition rate. However, at 
higher temperature, the regime shift to a diffusion controlled. But operation at higher 
temperature resulted in poor fiber quality and morphology which is elaborated in the 
latter part of this report. 

Least square analysis of the deposition rate gave the value of k to be 4900. Hence, 
the overall deposition rate is expressed as : 


k = 4900e kt 


C O.8 iO0.2 

BCl^TiCU 


f 


where k is the deposition rate in mm/hr, 

T is the temperature in Kelvin, 

Csch an< i C TiCl, are concentrations in mol/lit, and 
R is the gas constant =1.987 cal/mol K. 

This is the simple expression obtained for our operation. However, actual reaction 
kinetics is much more complicated due to many side reactions and products. 


7.1.3 Morphological Study 

Chemical vapor deposited TiB 2 exhibits various morphologies and microstructures de- 
pending on the temperature of the substrate and concentration of the reactants. It 
is advantageous to study and understand this phenomenon with respect to the final 
property of the fiber since the mechanical properties of the fiber is sensitive to the 
microstructure of the deposit. 

Temperature is one of the important parameter which determines the deposit char- 
acteristics. It is well established that at fixed concentration of reactants , increase in 
temperature, increases the size of the crystallite on the deposit. Vice versa also holds 
good [Besmann and Spear, 1975; Takahashi and Kamiya, 1974]. That means temper- 
ature and supersaturation play opposite to each other in determining the crystallinity. 
These effects are explained below. In Figure 12, the detailed surface morphology of the 
deposit at different temperature is shown. 
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The BCl 3 /TiCl 4 = 2 and H 2 /TiCl 4 = 22 was maintained throughout the opera- 
tion. At 1200°C, the coating (deposit) is smooth, microcrystalline without any nodular 
deposit. As the temperature increases, the surface becomes increasingly nodular. At 
1400°C, the surface is seen to be made up of nodules of over 50 /im in size which in turn 
made up of small grains of 2 to 5 fim size. After careful observation, it is predicted that 
initial deposit is microcrystalline and nodules appear to grow as a secondary features on 
it. At very high temperature (>1450°C), the deposit morphology changes from nodular 
to platelet kind of deposit which can be easily detached from the core [size around 30 
x 40 fim]. This may be due to the change in transfer mechanism [changed to mass 
transfer control]. At low temperatures (<1200°C), the deposit is amorphous in nature 
and hardly sticks to the substrate. 

The effect of concentration of reactants on the surface morphology is also consid- 
ered in our study. Figure 13 shows the effect of BClz/TiCl 4 ratio at optimum temper- 
ature (1250°C), flow rate (2 1/min) and H 2 /TiCl 4 ratio (25) on deposit morphology. It 
is observed that, increase in BCl 3 /TiCl 4 ratio increases the surface crystallinity leading 
to more nodular surface. At BCl 3 /TiCl 4 ratio greater than 2.5, nodules of approxi- 
mately 10/xm size are uniformly distributed on the surface of the fiber. However, at 
BCl 3 /TiCl 4 ratio of 2 or below the surface is smooth and perfectly microcrystalline. 

The effect of H 2 /TiCl 4 ratio on fiber morphology is presented in figure 14. The 
temperature and flow rate were maintained same as above, where as BClz/TiCl 4 ratio 
was kept at 2. At H 2 /TiCl 4 ratio of 7, the surface is very rough with a lot of small 
nodular growth, which are non-uniformly distributed on the surface. As the ratio in- 
creases, this phenomenon starts disappearing. At a ratio of 22, the surface is smoother 
with preferred crystallinity. Further increase in H 2 /TiCl 4 ration, no appreciable change 
in the surface features are observed. 

Even though we studied all possible parameters effect on morphological feature 
of the deposit, it is very difficult to locate the exact mechanism of growth. It may 
be influenced by impurities (side products), codeposited products, or any of the above 
parameter. However, optimum conditions for best quality microcrystalline fiber are at 
1250°C with H 2 /TiCl 4 and BClz/TiCl 4 ratios 25 and 2 respectively. 
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7.1.4 Thermal Stresses in the Deposition Process 

The fibers and their coatings are being exposed to three different stresses: thermal 
expansion mismatch, thermal stresses and blistering. These three stresses induce micro- 
cracks in the coating with the subsequent weakening of the final product. The question 
to answer is which one of these stresses crosses the resistance threshold for the cracks 
in the coating to appear. For most situations, a slow cooling process suppresses the 
thermal stress. On the other hand, blistering most likely will appear when coating and 
core interact to form a tertiary intermediate layer. The solution to thermal expansion 
mismatch is by incorporating intermediate layer which act as a ’’cushion” and absorbs 
stresses on the coating. 

The cooling rate df the deposited fiber plays also an important role on the ad- 
herence of the T 1 B 2 film on W core. If the fiber is brought to room temperature very 
fast, the coating chips off as shown in figure 15. This is due to the thermal stresses 
created while depositing. The slow cooling process after the reaction will suppress these 
thermal stresses leading to good adhesion. The parent fiber with the chip is shown in 
the figure. However, chipped piece showed remarkable good physical and mechanical 
properties. In our experiments, we limited the interaction between T 1 H 2 deposit and 
the W core to avoid blistering. SEM studies were conducted on a radial cross section 
of the fiber to determine possible interaction. 

There is no interaction and a sharp interface is clearly discernable. To further our 
studies, we mounted the fiber longitudinally. The special technique has to be used here 
to avoid the damage to the fiber and for accurate results. Fiber is mounted in an epoxy 
and polished first with alumina and subsequently with diamond wheel. These results 
axe shown in figure 16 where there is absolutely no interaction. The fiber is powdered 
and x-ray diffraction was carried out. However, there is no trace of tungsten borides as 
can be seen from x-ray pattern (figure 17). 

An intermediate layer of pyrolytic carbon was deposited to improve further the 
mechanical properties by eliminating blistering effect. Earlier analysis were used to 
detect the interaction between the core and coatings. Figure 18 does not indicate the 
presence of any diffusion layer or compound formation between the different layers of 
coating. This is in agreement with the results of Pierson et al., [1979]. 

Finally, the solution to mismatching thermal expansion coefficients has been to 
deposit an intermediate pyrolytic carbon layer; this ’’spongy” layer acts as a ’’cushion” 
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Figure 17: X-ray Diffraction Pattern of TiB 2 on W 


and "absorbs” stresses on the coating. However, the pyrolytic carbon layer may have 
a deleterious effect on the final strength. Figure 19 shows a layer of pyrolytic graphite 
deposited on W core at 1000°C. 

It can be seen that the external surface of graphite was striation corresponding 
to W core. However, final TiBi deposit shows the very smooth surface. That means 
the pyrolytic carbon (graphitic) layer provides a healing effect to the fiber stresses by 
smoothening the surface of the core. 

7.1.5 Residual Stress Model 

A theoretical model is used to estimate the residual stress in CVD of T 1 B 2 on W 
core. This provides information on the magnitude of residual stress and the region 
of maximum stress. The effect of various parameters can be studied and utilized in 
reducing the stresses. Since the thermal expansion mismatch plays a major role on 
residual stresses, only this is considered in our study. These stresses arise due to cooling 
the coating from deposition temperature. However, for high temperature application, 
for which the fibers are intended, these stresses are expected to diminish. 

A composite cylinder model is used [Brugger, 1971; Redston and Stanworth, 1945; 


37 





38 


ORIGINAL PAGil IZ 
OF POOR QUALITY 




39 


jfi.tr- 





Jain and Wetherhold, 1989], where coating and the core axe assumed to behave elas- 
tically in the temperature range considered. Due to anisotropy of the coating, it is 
modified to take into account the different properties of the coating along the different 
axes. The coordinate system used is shown in Figure 20. There are three components of 
stress in the cylindrical composite: axial stress c z , radial stress a r ,and circumferential 
or hoop stress cr fl . The constitutive equations and derivations are given in the Ph.D. 
thesis of Mr. Vikas Lakhotia, SUNY Buffalo (1990). 

The stress distribution in composite fiber is shown in Figure 21. The properties 
of the W were obtained from the literature [Shaffer, 1964]. The stress distribution is 
calculated for cooling a TiB 2 on W fiber from 1200°C (deposition temperature) to room 
temperature. Final diameter of the fiber assumed to be 200/im on 75/rm W core. It 
r an be seen that the stresses are compressive in the core and tensile in the coating. 
This is due to higher thermal expansion coefficient of TiB 2 compared to W. The same 
calculations, assuming the coating to be isotropic are given in Figure 22. 

There is not much difference except in the theta stress. In these calculations, it 
was assumed that the coating exhibited properties of polycrystalline TiB 2 . From the 
above figure, it is seen that the core is subjected to high compressive axial stress. In 
the coating, the radial stresses are compressive and hence not critical in determining 
the fracture point of the fiber. The theta and axial stresses are tensile and have a high 
positive value. The theta value is maximum at the interface and decreases rapidly with 
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radial distance. 

It is interesting to study the relationship between the core diameter to the final 
fiber diameter with respect to stresses. So we varied the p (fiber diameter to core 
diameter ratio) to calculate the generated stresses from the model. Since the maximum 
stress occurs at the coating interface, this is plotted in figure 23. It is seen that by 
increasing the p, decreased the stress and their levels off for the p value above 3. Above 
this ratio, stress does not change appreciably. 

The effect of deposition temperature on the residual stress is given in figure 24. 
Here the residual stress increases with temperature, as expected. This is obvious because 
higher deposition temperature means, greater the thermal expansion mismatch, so the 
fiber should be deposited in an optimum temperature to get microcrystalline deposit. 
The effect of internal stress on the stress field is also studied and is discussed in the 
later part of this report. 

7.1.6 Mechanical properties 

It is necessary to study the mechanical properties of the fiber to determine the optimum 
conditions for the deposition to synthesize a fiber. Measurement of the mechanical 
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properties of brittle material like TiB 2 fiber is a difficult task and often the technique 
used can affect the measured strength. All the tests were conducted on a research grade 
INSTRON tensile testing equipment. The technique used is the currently accepted 


standard for testing the brittle fibers. 

All the fibers studied fractured in a brittle manner. No separation of the coating 
from the core is observed. The typical fracture behavior is shown in figure 25. The crack 
directions on this photograph can be explained in terms of residual stress. This is shown 
in Figure 26. In the figure, external tensile stress is plotted against stress distribution. 
As the tensile stress increases, the theta and radial stresses remain unaffected, whereas 
axial stresses increase rapidly. In general, stress in any direction, causes the fracture in a 
direction perpendicular to it. Hence in the figure, the fracture is in a plane perpendicular 
to the fiber axis. Some cases, the cracks are in the radial direction due to reasonably 
high tensile stresses in theta direction. 

The effect of temperature on the strength of the fiber is shown in figure 27. The 
strengths measured is in the range of 30-50 KSI which is much below the hot pressed 
TiB 2 (100-150 KSI) and is pass through the maximum at around 1550K. This can be 
explained by considering both the residual stresses and surface morphology (microcrys- 
tallinity). In all the cases the deposition conditions were kept constant (2 min.). From 
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Figure 25: A Typical Brittle Fracture of TiB 2 fiber 
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Figure 26: Effect of External Tensile Stress on Stress Distribution 
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Figure 27: 

Strength Dependence 

on Deposition Temperature 


earlier figure, it is observed that residual stress decreases with coating thickness. This 
enables the coating to sustain greater external load prior to fracture. Here one should 
remember the fact that the coating thickness increases with temperature. However, 
with increase in temperature, the surface micro-crystallinity is lost. The surface be- 
comes more nodular as reported in our morphological studies, resulting in a poor fiber 
strength. These observations can also be explained by external stress effect plot (figure 
28). It shows that an external stress of 2.5 x 10 8 Pa («36 KSI), causes axial stress in 
the coating and the interface, of 16.5 x 10® Pa (~120 KSI). That means the coating is 
capable of carrying much larger load provided, the residual stresses are reduced. 

Various techniques were considered. Among them is the coating of intermediate 
spongy layer of pyrolytic graphitic carbon. The function of this layer is not clearly 
understood, but is believed to serve as a spongy layer capable of accommodating the 
stresses due to thermal expansion difference. The effect of pyrolytic graphitic carbon 
thickness is shown in figure 29. The tensile strength of the fiber due to the presence 
of the intermediate layer increases with the thickness. It is passing through maximum 
at around 3.5 \i m for 200 /im TiB 2 fiber. This gives a ratio of intermediate layer to 
final coating thickness ratio around 60. The modulus of the fiber is around 20-30 MSI, 
which is not affected by intermediate layer. 
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Figure 28: Effect of External Stress 






7.1.7 Texture of the Coating 


TiBi is a hexagonal crystal of AIB 2 type. The lattice parameters are, a = 3.028A, c - 
3 228A and -= 1.064A. It exhibits considerable anisotropy in thermal, mechanical and 
electrical properties along the two different axes, which is reported in literature [Mersol 
et al., 1968; Vahldick, 1967; Keihn and Keplin, 1967]. 

Many researchers observed preferential crystal orientation in TiB 2 deposit [Gan- 
non et al., 1963; Gebhardt and Cree, 1965; Besmann and Spear, 1975]. However, all 
the studies conducted are on flat substrate where C- axis of the crystal is parallel to 
the surface. This determination of orientation is fairly simple compared to on a fiber 
substrate. Following technique is used. 

The fibers were cut transversely, mounted in bakelite and polished by alumina 
and diamond. The coatings were examined metallographically. The samples were then 
etched with a mixture of nitric acid and hyrofluoric acid. The resultant photomicro- 
graphs taken were not able to conclusively prove the existence of orientation of the 

crystals. 

An X-ray diffraction technique described by Culhty was used [Cullity, 1978]. A 
Debye Scherrer pattern and a pinhole transmission pattern were produced using Ni 
filtered Cu-Ka radiation. From the diameter of the circle the angle of reflection (29) 
and the d-spacings can be found. Also the d-spacings for the hexagonal TiB 2 crystal 
are calculated and the results are given in Table 2. 

It can be seen that the [001] reflections are absent. This means that the [001] 
planes are perpendicular to the fiber axis. The Debye Scherrer is shown in figure 30. 

The results are confirmed by a pinhole transmission pattern, as shown in Figure 
31. The angle 8 is calculated from the ring diameter and 8 is measured from the normal 
as indicated in the figure. Using these two, the angle <j> is calculated by applying the 
expression, cos (f> = cos 8 cos 8. These axe given in 3. 

Table 4 gives the interplanar angles between the various planes. It can be seen by 
comparing the angles that [001] are the indices of the fiber axis. This means that the 
c-axis of the TiB 2 crystal is parallel to the W-fiber surface. The second ring suggests a 
second, much less pronounced [Oil] texture. 

TiC and B 4 C were also prepared by CVD. These fibers were obtained by decom- 
position of the corresponding halides with a hydrogen and carbon source. 
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Table 2: Debye Scherrer Pattern 
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20.0 
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1.020 

20.2 

1.015 

10.3 

0.992 

21.1 
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0.8774 

21.2 
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Table 3: Pinhole Transmission Results 


1391 

e 

6 


Qj 

17.1 

90 



17.1 

63 

64 

^^9 

17.1 

113 

112 

10.1 

22.3 

46 



48 









Table 4: Interplanar Angles 



8 Boron carbide and titanium carbide fibers 

[See Figure 32 and 33]. 

The CVD process of B A C is carried out on a carbon core by reducing BCl 3 and 
CHi/CCU with hydrogen in the temperature range 1000-1500°C. The concentration of 
BCh, hydrocarbon (CH t / CCl 4 ) and hydrogen gases was varied in our experiments. 
The nucleation of excess of carbon was observed when CH$ was used as a precursor. 
This nucleation of carbon within the deposit is catalyzed by BCI 3 as reported earlier 
[Brennfleck et al, 1984]. This process leads to local separation of B 4 C and carbon 
deposition. This intergranular carbon deposition exhibits a very low Young’s modulus 
and enables elastic strains of the whole deposit without plastic deformation. The prob- 
lem of low modulus was partially overcome by using CCl 4 as the carbon source. The 
optimum reactant concentration expressed in terms of Cl to H ratio was found to be 
1:5 at 1300°C. An excess of hydrogen increases the conversion per pass of BCl 3 and 
reaches the maximum at around 0.002 m 3 / sec (total gas flow rate). Our best deposit 
contained less than 2% of free carbon. Higher temperature gave many side products 
[BHCl 2 , BCl 2 , and BCl as reported by Naslain et al., 1979]. The results are given in 
Table 5. Our fiber has tensile strength of approximately 300-400 KSI with a modus of 
approximately 35 MSI. The fiber contained less than 2% of free carbon. 

Titanium carbide deposition was carried out by reducing TiCl 4 + CCl 4 with hy- 
drogen in the temperature range 1300-1500°C. The optimum conditions were, chlorine 
to hydrogen ratio 1:7 with a flow rate of 0.0015 m 3 / sec (at 1400°C) and the residence 
time in the reactor is 3 min.' Here we are faced with a severe problem of heterogeneous 
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Table 5: Properties of fabricated fibers 
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Temp, of 
deposition 

(°C) 
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xlO 6 °C~ 1 
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Oxidation 

resistance 
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carbon 

| 

3180 

1 
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<4% 


2450 


5.0 


25-35 

Poor 

<2% 

1 

2980 


7.9 


20-30 

Moderate 


■Efl 

2000 




30-50 

Poor 



nucleation and simultaneous deposition of carbon and TiC. In the CH A /TiCl A g as ma- 
ture, the TiC nucleation is never inhibited if an excess of TiCl A i s used. Carbon core 
strength (fiber strength) decreases by coating. This decrease is due to the diffusion of 
carbon from the core substrate to the surface because of the inhibited deposition of 
carbon from the gaseous phase. Initial analysis showed that the best result for the de- 
position process can be obtained with a carbon to titanium ratio ^3. Our fibers showed 
a tensile strength of approximately 200-300 KSI with a modulus of 35 MSI. Free carbon 
in the mixture is less than 5%. Results are tabulated in Table 5. 


9 Chromium Diboride Fibers (CV.B 2 ) 

9.1 Introduction 

Chromium diboride should be very effective in resolving or reducing severe mechanical 
and high temperature corrosion problems, especially in severe environments. Further- 
more, its ductility, density, good physio- chemical compatibility with metals, high ex- 
pansion, all make it one of the best candidates among the refractory materials for the 
metal-matrix composites in the space industry. 

The CVD of chromium is well known and the structure-property-process rela- 
tionship have been extensively studied. [Mazille, 1980; Wakefield, 1969, Hanni and 
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Hintermann, 1977] However, very few or no investigations have been made using a 
CtB 2 deposition either on a fiber substrate or any other shaped substrate. The most 
important aspect in CVD of CtB 2 is obtaining the chromium species in the gaseous 
phase at a relatively low temperature and with ease. Different precursors have been the 
subject of many chromium deposition studies in the past, but only few of them have 
met the above mentioned requirement. These previously investigated precursors are 
CrCl 2 [Hanni and Hintermann, 1977], Cr 2 O z [Roos, 1959], Cr(CO) G [Owen and Web- 
ber, 1948], Organo-chromium compounds [Wakefield, 1969; Oxley, 1962] and chromium 
sulphides. Except for Cr(CO ) 6 and the organo-chromium compounds, the rest of the 
above mentioned precursors have either a high melting and boiling point or handling 
difficulty. The carbon content of the lower melting compounds [Cr(CO ) 6 , (organo- 
chromium compounds)]' lead to the formation of carbides and oxides which are stable 
compounds at the temperature of deposition. Considering the availability, ease of han- 
dling and the high vapor pressure at relatively low temperature [M.pt -96.5°C, B.pt 
116. 7°C] of chromyl chloride [Cr0 2 Cl 2 }, it was decided upon as the chromium con- 
taining precursor in this study. Moreover, and as far as we know, no one has been 
specifically devoted to the in-depth CVD of CrB 2 using this precursor or any other 
precursor. 

In the present work we have used a thermodynamic approach based on the min- 
imization of the total Gibbs free energy of the system, to study the influence of the 
CVD parameters on the composition of the deposit obtained at equilibrium, from the 
Cr0 2 Cl 2 + H 2 + BCl z reactants system. The emphasis of this study was to obtain 
pure CtB 2 as a condensed phase. Several experiments were performed to study the 
morphology of the deposit. 

Modelling the kinetics of any chemical process is extremely difficult. This is es- 
pecially true for deposition processes. In these processes a large number of reactions 
take place simultaneously, being either gas phase or surface reactions. A classical way 
around this has been to identify the rate limiting reaction step and then modelling this 
step accurately. As the rate controlling step is dependent on both the reactant flowrate 
as well as the reaction temperature, this approach, although sufficient in most cases, 
cannot be extrapolated for use in scaled-up systems. 

The modelling of CVD-kinetics with simultaneous reactions has been done[Coltrin 
et al., 1986]. Although encouraging results were obtained by this method, a parameter 
was always used to adjust the calculated deposition rates to the experimentally observed 
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rates. This is done by adjusting the apparent activation energy of the surface reactions. 
A detailed description of the flow profiles in the reactor is needed for the theoretical 
approach. As flow profiles tend to be three dimensional through a horizontal reactor, 
even with a symmetrical set up, this translates into a very complex problem. [Scholtz 
1991]. Experimental work that had been carried out had used a Freundlich isotherm to 
describe the deposition reaction. 

0 = C * P* (3) 


6 Surface coverage 
C Rate constant 
P Adsorbing gas Pressure 
n Real >1.0 

Results obtained from such a study resulted in the curve fitting of the experimen- 
tal data to a predicted rate equation and did not describe the actual mechanism of 
deposition. The importance of the different precursors can be identified with such a 
kinetic expression, but any other detail of the reaction mechanism cannot be deduced 
from the fitted rate equation. It is important to note that the Freundlich isotherm can 
be used to approximate most compiled reaction mechanisms. Keeping all these factors 
in mind a different approach was decided upon to model the reaction kinetics. 


9.2 Thermodynamic Analysis for CtB 2 from CrO^Cli + H 2 + 
BClz CVD system 

A thermodynamic approach has been used to study CVD of the CrOzClz- # 2 “ BCI 3 
system. But, in principle, this may not necessarily give information on the CVD- 
mechanism or on the deposition characteristics. It can however be used to get a broad 
overview of the system, assuming that states close to equilibrium are reached in the 
vicinity of the substrate. In general these conditions may not be reached in each and 
every case. 

The thermodynamic approach which is used here and its application to CVD is 
based on the minimization of the total Gibbs free energy of the system as described 
elsewhere. [Colmet et al, 1982; Christin et al, 1979; Zhao et al, 1990]. This method 
gives the solid, liquid and vapor phase composition at equilibrium for a known initial 
composition, temperature and pressure. Here one should assume that the deposit is in 
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equilibrium with different phases and that all the chemical species able to be present 
are known. 

The CVD oiCrB -2 from a Cr0 2 Cl 2 - H 2 - BCl% mixture is based on the reduction of 
Cr0 2 Cl 2 and BCl 2 by H 2 at the reaction temperature. The actual reaction mechanism 
may be very complex whose result strongly depends on the deposition parameters. 

In the present investigation 53 different possible chemical species were considered: 
eight in condensed phase (CrB, CrB 2 , CrCl 2 , CrCl 2 , Ct 2 0 2 , Cr) six in the liquid phase 
and the rest of them in vapor phase. The values of their thermodynamic properties 
(Gj.- #298 ) were taken from tbe JANAF table [1970] and Handbook of Thermophysical 
properties [Barin et al, 1977]. The concentration of all the species formed at equilibrium 
were investigated as a function of temperature (800 K < T < 1500 K), total pressure 
(0.5 atm. < P < 2.0 atm.) and initial composition of the vapor phase defined by the 
a= H 2 /Cr0 2 Cl 2 (1 < a < 200) and 0 = BCh/Cr0 2 Cl 2 (1 < 0 < 10) ratios. 

At a given pressure (1 atm), the nature of the chemical species depends on the 
temperature a well as the initial reactant concentration of the vapor phase. In order to 
identify an optimum operating region we varied the a ratio with temperature, keeping 
0 constant at 5. The value for 0 = 5 is justified in our later observations. Figure 34 
gives the product distribution at three different temperatures with a. 

At high a values, the thermodynamic yield of CrB 2 remains constant irrespective 
of temperature. However, at low a values, CtB 2 increases with temperature before 
stabilizing. The rate of formation of CtB 2 is highest at low temperature and a values. 
One of the important condensed species which has an authoritative effect on the material 
property is free elemental boron. It is observed at high temperature (above 1300 K) for 
all values of a and increases with temperature. 

A study in the variation of the product distribution with a variation in the BCl 2 
initial concentration is very important as it will definitely give us some insight in the 
B to Cr ratio necessary to produce CtB 2 . All calculations were carried out at a= 100. 
From figure 35 it can be seen that only CtB 2 , CrB, B and CrCl 2 are thermodynamically 
stable condensed species. 

It is interesting to observe the parabolic profile for CrB 2 , which is passing through 
a maximum at a different level, with respect to temperature and the 0 ratio. Free 
elemental boron is prominent at a high 0 ratio and a high temperature. At lower 
temperature (T = 1100 K ) , ' i t is present in negligible amount. Chromium dichloride is 
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observed at low temperature (< 1100 K) and low f3 ratios. This completely disappears 
at high values of /?. 

Keeping all this in mind one can conclude that the mechanism of deposition is very 
complex depending on temperature, a ratios as well as ratios. At low temperature 
and low values of a and /3, the chromium precursor goes through the CrCl 2 phase. This 
shows that initial mechanism of formation of chromium is possible from Cr 2 0 2 Cl 2 via 
CrCl 2 reduction. As the temperature increases and at high ratios of a this reduction 
is total. With an increase in operational values, the mechanism seems to change into a 
replacement reaction mechanism where BCI 3 reacts with chromyl chloride directly and 
the boron atom displaces the Cr atom from the chromyl chloride molecule. This leads 
to the formation of mostly of B 2 0 2 and B 2 0 2 C1 2 at high values of with no Cr0 2 
formation. From the overall view one can suggest that the majority of the CrB 2 is 
formed due to a reduction reaction. This leads us to have the following values for our 
experimental operation: a should be 90 < a < 110, (3 should be 4.5 < (3 < 6.5 and the 
temperature should be between 1100 and 1500 K to get a maximum amount of CtB 2 
without any impurities. To get an overall pictorial view, a three dimensional plot was 
drawn [figure 36]. 

The effect of pressure on our product distribution was studied by varying the 
pressure and the a ratio, while keeping both the (3 ratio and temperature constant at 5 
and 1300 K respectively [Figure 37]. Three different pressures (0.5 atm., 1.0 atm. and 
2.0 atm) were used and their product distributions were compared with each other. A 
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Figure 37: Effect of Pressure on Product Distribution 

lower pressure increased the formation of all the condensed phases. (CrB, CrB 2 , B, 
B 2 0 3 (l)) This resulted in the formation of CrB 2 at a lower temperature as well as at a 
lower a ratio. Unfortunately the formation of B also increased malting lower operating 
pressures not more attractive. Although the CrB 2 formed at lower temperature and 
a ratio, the ultimate conversion of boron into CtB 2 stayed the same irrespective of 
pressure. 

9.2.1 Results and Discussion 

By looking at the phase diagram of chromium and boron, it is apparent that in order 
to get CrB 2 as the most stable condensed phase we should keep the boron to chromium 
ratio greater than two. [figure 38]. 

We divided the borides broadly into two categories, lower borides (CrB) and higher 
borides ( CtB 2 ). From the experimental results we can however detect many multiple 
borides such as CrB 4 , Cr 2 B, Cr 5 B 3 , etc. by X-ray diffraction spectrograph. How 
critical these multiple borides are in application to composites is yet to be seen. The 
typical analysis of the results are discussed below. 

The experiments were first carried out at thermodynamical optimum conditions 
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to get exclusively higher borides as the condensed phase. As far as possible we tried to 
avoid getting free boron in the deposit. The conditions that were used here were a= 
110 and j3= 5. Since microcrystallinity is also an important parameter for obtaining 
a better quality fiber, we investigated the deposit structure also with respect to the 
substrate temperature. Figure 39 gives the SEM photographs at a = 110, (5— 5 and 
various temperatures. At a temperature of 1230 K the deposit is more amorphous and 
loosely bound to the fiber. 
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With a slight increase in the substrate temperature, the deposit became very 
smooth and microcrystalline which is the desirable quality for a good quality fiber. 
This was observed around 1273 K. A further increase in the temperature leads to an 
increase in crystal size and irregular deposit. One can see the irregularity right from 
the initial point of deposition. This may be due to the impurity effect which can act 
as a ” catalyst” leading to a non homogeneous growth. Again at high temperature, the 
mass transfer effect should play an important role as the reaction rate is no longer the 
controlling deposition parameter. Whatever the reason may be, the fiber quality was 
very poor. Several runs were taken varying the a and (3 ratios with temperature and 
observing their effect on the microstructure a well as the composition of the deposit. 
It is observed that a higher (3 ratio results in the incorporation of more free boron in 
the deposit and a lower' a as well as (3 ratio suppresses the formation of free elemental 
boron but results in the formation of multiple borides ranging from CrB, CrB 4 , CtB 2 
to CV 5 .B 3 , At high temperatures (1553 K) a small amount of Cr^O^ was also observed 
in all the cases which contradicts our thermodynamic predictions. 

Experiments were also carried out to see how the precursors CrClz and Cr(CO ) q 
behaved with respect to chromyl chloride. The difficulty of obtaining chromium trichlo- 
ride in the gets phase, condensation in the feeding lines and uneven deposit limited its 
analysis. The chromyl hexacarbonyl deposition results at 1350 K are shown in figure 
40. This shows an uneven deposit which can be explained because of the large number 
of impurities in the deposit. From the x-ray diffraction pattern it was observed that 
mainly B 2 hC y Cr0 2) C r r 2 3 C r 6 , Cr and a small amount of CtB 2 was found in the deposit. 


Finally a strong penetration of deposit into the tungsten substrate wa£ observed. 
The quality of the fiber is reduced as the core is weakened by the infiltration reaction. To 
counter this problem some fibers with intermediate carbon layers were deposited. From 
the results it seems that this gave both an improvement in the quality and appearance 
of the fiber. This can be seen in figure 41. The fibers from these runs showed the tensile 
strength up to 167 KSI. The strength analysis was carried out by NASA scientist. 


9.3 Kinetic Study of the CrB<i Deposition Reaction 

As a first step in the prediction of a kinetic mechanism for the deposition reaction, a 
set of elementary reactions were predicted. This system of reactions was derived from 
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Figure 41: Fiber with Intermediate layer of Carbon 
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the previously found thermodynamically most stable species and taking into account 


the bi-molecular bond energies hetween the different atoms. 

H 2 ^ 2 H 

(4) 

BCl 3 + H 

^ BCl 2 + HCl 

(5) 

BCl 2 + H 

^ BCl + HCl 

(6) 

BCl + H 

^ B + HCl 

(7) 

Cr0 2 Cl 2 -b BCl 

^ BOCl + CrOCl 2 

(8) 

CrOCl 2 + BCl 

^ CrCl 2 + BOCl 

(9) 

' ZBOCl 

B 2 0 3 C1 3 

(10) 

CrCl 2 + H 

^ CrCl+HCl 

(11) 

CrCl + H + BCl 

^ CrBCl + HCl 

(12) 

CrBCl + H + BCl 

^ CrB 2 + HCl 

(13) 


Three different reaction models were proposed and their agreement with experimental 
results is tested. 

It is well known that the deposition reaction is a combination of both gas and 
surface reactions. The aim of this study was to identify the stage at which adsorption 
onto the surface, to form the deposit, occurs. It was observed thermodynamically that 
the chromium containing precursor Cr0 2 Cl 2 , in the presence of both H 2 and BCl 3 , 
is reduced to the previously identified compound CrCl 2 at elevated temperatures. It 
was found that 250°C was achieved in the preheater which was enough to complete 
the reduction reaction described above. It was thus decided to use the chromium 
intermediate CrCl 2 , generated in the preheater, as the starting chromium compound 
in the modelling of the deposition reaction. 

9.3.1 Results and Discussion 

All the experiments were done in the kinetically controlled regime. The first approach 
to test the model was by assuming a pure gas phase reaction and resultant CrB 2 solid 
deposit (adsorbs) on to the substrate. Taking all the gas reactions at equilibrium, a 
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simplified model was derived. The rate equation is written as 

Rate = ( 14 ) 

p 3 p 14 
r B 7 O z r HCl 


It is apparent that both the exit and entrance gases should be monitored to inves- 
tigate this model. However, our measured rate showed a very strong surface reaction 
dependence. Figure 42 shows the change in surface area of the substrate effect on 
the deposition rate. From the graph, it is evident that initial deposition is depend on 
surface area of the substrate. After certain time, both treated and untreated samples 
(having surface area 0.09 m 2 / gm, 0.05 m 2 / gm resp) showed the same rate. Hence our 
first model was not valid. 


The two different adsorption model were assumed to test our results: 2nd model, 
adsorption of the CrCl 2 onto the substrate was assumed with subsequent reduction by 
BClz and H 2 . Adopting the Langmuir isotherm, the rate equation is expressed as 


Rate = 


1 + KlPcrClj 


(15) 


Fined model (3rd) assumes that both BCl z and .CrCl 2 adsorbs on the surface and 
reduced in the presence of H 2 to deposit CrB 2 . This can be schematically represented 
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Figure 43: Schematic Representing the Model 


in the sketch (figure 43). The elementary reactions are: 


CrCh + S’ CrCl' 

(16) 

BCl 3 + S' ^ BCl 3 

(17) 

2 BCi; + 4 H 2 + CrCl 2 ^ CrB; + 8 HCl 

(18) 

CtB * ^ CrB 2 + S' 

(19) 


By assuming Langmuir isotherm, and both BCI 3 and CrCl 2 adsorbed species were 
equally mobile on the substrate, the rate expression is given as 

„ kj K 1 KlPc 1 c h no^ 

“ e (1 + K 1 P c ,ch + K,P Bch )' 

Fitting the experimental data to the proposed models created some problems. A 
non-linear least squares fit is used to find the constants that best approximated the 
experimental results. The following results were obtained. For the second proposed 
model (Single adsorption of the CrCli onto the substrate), no good fit for the experi- 
mental results could be found. Fitting the experimental results to the third proposed 
model gave much better match. The maximum deviation between the experimental and 
proposed model was 8%. The agreement between the experiments and the proposed 
model was plotted on a 3-D graph, to show the variation between the experimentally 
found deposition rates and the calculated deposition rate. [Figure 44). 

A higher resolution plape was generated in the three-Dimensional domain so that 
the variation of the deposition rate with a variation of the model variables could be 
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observed. At the same time the best fit for the second proposed model was also plotted 
and the difference in the generated graphs is apparent. 

To get a feeling of how good the agreement between the experimental and model 
results, a comparison between the obtained rates were drawn. 

Although the agreement is not excellent, this definitely shows that an adsorption 
model in which both the BCl$ and the GtGIi species adsorb on the surface describes 
the actual mechanism of the deposition reaction. Using a more refined model, where the 
surface reaction occurs in a few steps with the slowest step determining the deposition 
rate, should give better agreement between the experimental and theoretical models. A 
much improved experimental set-up would be needed to be able to do a more detailed 
kinetic study. 

i 

9.3.2 Implications of the Calculated Rate Constants 

Agreement between the experimental and the predicted dubbel adsorption model was 
quite good and the largest error was of the order 8%. The rate constants were de- 
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termined by using a least squares fit for the non-linear rate equation. The reaction 
coefficients were found to be the following values. 

ki = 113.0 (fim. min" 1 . atm ~ A ) 

K\ = 8.0*10 6 

K 2 = 4-0* Z0 6 

These values were found to best represent the experimental results and their rel- 
ative values were found from the least squares fit used to approximate the non-linear 
predicted deposition rate. Evaluating the order of magnitude of the constants, might 
give insight into adsorption efficiency as well as the reaction mechanism characteristics. 

10 CONCLUSIONS 

The CVD technique was successfully employed to synthesize TiB 2 , CrB 2 , B, TiC and 
B 4 C fibers on tungsten and carbon substrate. These fibers showed good physical and 
mechanical properties for the possible use in metal-matrix composite. The main pro- 
cess variables are temperature, compositions of the reactants and total flow rate. The 
thermodynamic studies on TiB 2 and CtB 2 , were used to predict a window for pure 
product. The rate of deposition of TiB 2 was expressed as: 

k = moe^cr Ch c° T ? cli ( 21 ) 

whereas CtB 2 deposition rate was best expressed by using an adsorption model as- 
suming both the reactant species ( CrCl 2 and BCl 3 ) adsorb on to the surface and get 
reduced with hydrogen to produce CtB 2 . 

morphological studies on TiB 2 and CtB 2) showed the optimum deposition condi- 
tion. The metallographic and X-ray diffraction studies on TiB 2 gave the TiB 2 crystal 
orientation parallel to the C-axis. The stress distribution due to thermal expression 
mismatch was studied on TiB 2 deposit and is partially neutralized by using interme- 
diate spongy coating of carbon. The similar experimentation on CtB 2 also produced 
better results. The tensile strength with TiB 2 reached up to 220 KSI where as with 
CrB 2 is 167KSI. Exploratory work with TiC, B^C and B fibers produced the strength 
result that are comparable to the reported values in the literature. Fibers showed good 
potential for the use in metal-matrix composite. 
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11 OVERALL SUCCESS OF THE PROJECT 


So far we have successfully developed the CVD process for synthesis of TiB 2 , TiC, 
B^C , CtB 2 and B fibers and the regimes of deposition, quality of synthesized fibers and 
qualitative characteristics of the process has been pin pointed. The TiB 2 fibers showed 
extremely good strength and continuous process fibers (several meters) has been sent 
to NASA for testing and analysis. With proper encouragement, the process is ready to 
work continuously. The TiC, B 4 C , B fibers matched the properties of the best known 
reported values. Our exploratory work on CtB 2 gave unexpectedly excellent results. 
Recent results from NASA showed for non-primary carbon coated fiber, as high as 170 
KSI. This process needs more work to go for full scale pilot plant. Overall, this project 
is fairly successful in developing the fibers with a potential for larger scale production. 
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13.1 Titanium Diboride Fibers 




Titanium diboride 

Formula weight : 

Melting point : 

Density : 

Chemical 

Theoretical analysis 

Reactivity : 


TiB 2 

69.54 g/mol 
2980 °C 
4.52 g/cc 

31.12 % boron 
68.88 % titanium 

with HC1, HF. 


1. Does not react 

2. Reacts with hot H 2 S0 4 . 


Reference 

82 


e &- 


3. Soluble in HN 03 -H 2 0 2 . 
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4. Decomposed by fused alkali, carbonates and bisulfates. 

5. Reacts violently with Na 2 0 2 and brown lead oxide. 

6. Air oxidation becomes noticeable at 700 "C. 
Electrical 


83 

84 


Resistivity : 9.1 fiD- cm. at 300 °C. 

The temperature coefficient of electroresistance is 0.84 x 10~ 3 deg -1 


Mechanical 


85 


Knoop Hardness 
Transverse Rupture Strength 
Fracture Toughness, K ic 
Youngs Modulus 
Poissons Ratio 

Thermal 

Conductivity 

Specific Heat 


15-36 GPa 
700-1000 MPa 
6-8 MPa m 1 / 2 
514-574 GPa 
0.09-0.28 

65-120 W/m K at 300 K 
54-122 W/m K at 2300 K 
632 J (kg K)’ 1 at 300 K 
1155 J (kg K)" 1 at 1400 K 


% 
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Expansion 

// a-axis : AL/L 0 = - 0.114 + 3.092 x 10" 4 T + 2.806 x 10" 7 T 2 - 3.996 x 10" 11 T 3 % 
// c-axis : AL/L 0 = - 0.197 + 5.959 x 10“ 4 T - 2.770 x 10~ 7 T 2 - 4.369 x 10 -11 T 3 % 
polycrystalline : AL/L C = - 0.142 + 4.054 x 10 4 T + 2.786 x 10 T 

- 4.097 x 10~ n T 3 % 


t 



RESULTS AND DISCUSSION 


BCl 3 /TiCl 4 =2.0 



Temperature^ 14 00 K 

Temperature^ 1500 K 

- — Temperature=1600 K 


j l— 

0.0 15.0 


20.0 25.0 30.0 

H 2 /TiCl 4 

t V 


Effect of H 2 /TiCL, ratio on Deposition Efficiency 



DESIGN AND DEVELOPMENT OF CVD UNIT 


TABLE 3.1 

Physical Properties of Tungsten (W ) 61 


Melting Point 

3417 °C 

Boiling Point 

5930 °G 

Vapor Pressure (mm Hg) 

1527 °6 

1.93 x 10- 15 

Specific Heat at 20 °C ( cal /gm- atom) 

6.25 

Thermal Conductivity (cal/sq cm/cm/sec^C) 

to 

o 

0 

O 

0.310 

CO 

to 

0 

Q 

0.275 

1127 °C 

0.268 

1527 °C 

0.253 

Density (g/cc) 

19.17 

Coefficient of Linear Expansion 

mean value 0 - 500 °C 

4.98 x 10" 6 ( 0 C)- 1 

Elastic Properties (at 25 °C) 

Young’s Modulus 

41 x 10 3 kg/mm 2 

Shear Modulus 

16 x 10 3 kg/mm 2 

Poissons Ratio 

0.27 

Electrical Resistivity (/ifi-cm) 

to 

0 

O 

5.89 

o 

o 

0 

o 

22.43 

1100 °c 

34.65 

1500 °C 

49.66 
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Metallograph of CVD TiB 2 
















A. Ti + BC1 3 + H 2 at 1400°C 

B. Ti fiber + boron powder + hydrogen carrier gas at 1350°C 

C. Ti fiber + BC1 3 + H 2 + Ar at 1200°C 







L. LOOO°C 
3. I L 50 ° C 


L 2 2 0 ° C 


03&1NA- ’.'AGE 'S 
OF POOR QUALITY 

H : C 1 ratio • 7 

B : T i ratio •- 




KV X 450 1 0 0 U 106 51587 UBE 



dA 



I • >'i> j. 



(O) i > cc^e 
C b j Wo Cote 
( C y Ti Coie. 



°f> { 


• ;n 


COncb I *' ciXQ 


£,L>1Cxin 

ORIGINAL PAGE e'S 

OF POOR QUALITY 


C' cr* W 







30KV X 120 1 00U 403 70589 UBE 


9. UOO°C LO. U50°C 

11. >U50°C 

ORIGINAL PAGE IS 
OF POOR QUALITY 


30KV X340 100U 102 63087 UBE 







5. L250°C 

7. L 3 3 0 ° C 


6. 

8. 


1 2 50 ° C 

1 3 50 ° C 


OwGiNAL PAGc is 
OF POOR QUALITY 


tr: 






ORIGINAL PAQF i S 

F. L 2 50 ° C 0F POOR Q UALrr y 




J. 


greater than 1400°C 


Of F ’ M r#, i .rry 


I. 


1400 °C 




ORIGINAL PAGE IS 
OF POOR QUALITY 


/tv 


30KV X220 1 00U 118 60288 UBE 














WITHOUT INTERMEDIATE LAYER 


ORIGINAL PAGE IS 
OF POGJ? O’ '.*» rry 
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39KV X57^ "^05^^0^ 70589 UBE 
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OF POOR QUALITY 


TiC FIBERS ON CARBON CORE 

(a) TiCL^ deposition with CCl^ at I200°C with C:Ti = 3; H:CL = 7. 

(b) insitu conversion Jt L300°C 

(c) TiCL. + CCL, deposition at L300°C with C:Ti = 3; H:CL = 7.5 


I Cr‘ 





Til?2 with graded deposition of carbon 



30KV X370 
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c 

c 

AUTHOR 

: VIKAS 

LAKHOTIA 


c 

DATE 

: OCT . , 

1989 



THIS IS A PROGRAM USED TO CALCULATE THE THERMODYNAMIC C 
CONSTANTS al . . a7 FOR THERMODYNAMIC EQUILIBRIUM CALCULATION . C 
THESE CONSTANTS ARE USED IN THE PROGRAM NASAF TO OBTAIN C 
EQUILIBRIUM COMPOSITION OF REACTING MIXTURES. THIS PROGRAM C 
USES SEVERAL SUBROUTINES AVAILABLE IN IMSL LIBRARY C 
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PROGRAM LSQ 

IMPLICIT DOUBLE PRECISION (A-H, O-Z) 

INTEGER LDA, NOBS, NPARM, NEQNS , LDB, NB, NCA, NRA.LDC, 

* NCZ,NCC,NRZ,NRC 

PARAMETER (NOBS=8, NPARM=5 , NEQN=3 , LDA=NOBS*NEQN , 

> 

* LDB=NPARM , NB=NPARM , NCA=NPARM , NRA=NOBS*NEQN , 

* LDZ=NOBS*NEQN , LDC=NPARM , NRZ=NOBS *NEQN , NCC=1 , 

* NCZ = 1 , NRC =NP ARM , IPATH= 1 ) 

REAL *8 A(LDA,NCA) ,B(LDB,NB) ,Y(LDA) ,Z(LDA,NCZ) , 

* C(NPARM.NCC) ,D (NPARM) ,F (NPARM) ,X(NOBS) ,DHF,R, 

* Q(NOBS) ,S298,SQ(N0BS) .SUM.SUM1 
CHARACTER*32 FILNAM 

EXTERNAL DMXTXF , DMXTYF , DLSARG , DWRRRN 


write(* ,*) ’please enter the filename from which data is 
* to be read’ 

read(* , 100) filnam 
100 format (a32) 

open(l3 ,f ile=f ilnam, status=' old’ ) 
open(l4,file=’lsql .dat’ , status= ’unknown’ ) 
ne=0 


R=8 . 314 
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R is the gas constant 
read(13,*) dhf 
read(l3,*) S298 

dhf is the heat of formation of a substance at 298 K. 

S298 is the entropy at 298 K. 

This loop reads the value of specific heat, entropy and 

> 

enthalpy at any temperature, t, and arranges them in a 
column. It is used to calculate the five constants al . .a5. 
It calculates the elements of the matrix on 
the right hand side which when multiplied with the 
constants give the specific heat, entropy or enthalpy. 

do 56 i=l,nobs 
read(13 , *) x(i) 
t=x(i) 

read(l3 , *) y(ne+i) 

Here y(ne +i) is the specific heat. 

y (ne+i) =y (ne+i) /R 

A(ne+i,l)=l. 

A(ne+i , 2) =t 
A (ne+i , 3) =t**2 
A (ne+i , 4)=t.*'*3 
A (ne+i , 5)=t**4 
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ne=ne+l 

read(13,*) y(ne+i) 

Here y(ne+i) is the entropy. 

aa=y (ne+i) 

sq(i)=aa 

y (ne+i) = (y (ne+i) -S298) /R 
A (ne+i , 1) =log(t/298 . ) 

A (ne+i , 2) =t-298 . 
A(ne+i,3)=(t**2-298.**2)/2. 

A (ne+i , 4) = (t **3-298 . **3) / 3 . 

A (ne+i ,5)=(t**4-298.**4)/4. 
ne=ne+l 

read(l3,*) y(ne+i) 

Here y(ne+i) is the enthalpy 

aa=y (ne+i) 

q(i)=aa 

y (ne+i ) =y (ne+i) *1000. /R 
A (ne+i , 1 ) =t - 298 . 

A (ne+i , 2) = (t **2-298 . **2) /2 . 

A (ne+i , 3) = (t**3-298 . **3) /3 . 

A (ne+i ,4)=(t**4-298 . **4) /4. 

A (ne+i , 5) = (t'**5-298 . **5) /5 . 
do 57 i=l ,nobs*neqn 
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1 


aa=y(i) 

57 z(i,l)=aa 

C The following loop makes all the elements of the matrix of the 

C same order. 

do 61 i=l,lda 
A(i , 1) =A(i , 1) 

A(i,2)=A(i,2)*1.0E-03 
A(i,3)=A(i,3)*l. 0E-06 
A(i ,4) =A(i ,4)*1 .0E-09 
61 A(i,5)=A(i,5)*i.0E-12 

C The subroutine DMXTXF multiplies the transpose of the matrix 

C A with itself. The end result is a square matrix. 

CALL DMXTXF (NRA.NCA, A, LDA ,NB ,B ,LDB) 

C All the dwrrrn subroutines are used to display the results 

C of the calculations on screen and can be omitted without 

C any difference to the results. 

C CALL DWRRRN ( ’B=trans (A) *A’ ,NB,NB,B,LDB,0) 

C The subroutine DMXTYF multiplies the transpose of A with Z. 

C This is done to maintain equality of the two sides. 

CALL DMXTYF (NRA.NCA, A, LDA, NRZ,NCZ.Z,LDZ,NRC,NCC,C,LDC) 

C CALL DWRRRN ( ’ C=trans (A) *B ’ , NRC , NCC , C , LDC , 0) 

do 58 i=l,np‘arm 
aa=c(i , 1) 
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d(i)=aa 

This subroutine solves a system of linear equations. 

It brings back the value of the constants al..a5. 

CALL DLSARG(NPARM,B,LDB,D,1,F) 

CALL DWRBBH(»FM.HPARM,F,i.O) 

The following loop makes the elements of the correct magnitude 
and corrects for the changing of order of magnitude earlier. 
al=f (1) 

a2=f (2)*1.0E-03 
a3=f (3) *1 . 0E-06 
a4=f (4) *1 .0E-09 
a5=f (5) *1 .OE-12 

The folowing lines calculates the value of a6 using the 
enthalpy values. 
sum=0 . 

do 92 i=l .nobs 
t-x(i) 

q(i)=(q(i)+dhf ) *1000 . /R-al*t-a2/2 . *t**2-a3/3 . *t**3-a4/4 . 

* *t**4-a5/5 . *t**5 

sum= sum+q ( i ) 
a6= sum/ nobs 

The following lines calculate the value of a7 using the 


entropy values . 
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suml=0. 

do 93 i=l,nobs 
t=x(i) 

sq(i)=sq(i) /R-al*log(t)“a2*t-a3/2 . *t**2-a4/3 . *t**3- 
* a5/4 . *t**4 

93 suml=sumlH-sq(i) 

t 

a7=suml/nobs 

write (14 ,41) al ,a2,a3,a4,a5 
write (14, 42) a6 ,a7 

41 format (5el5 . 8) 

42 format (2el5 . 8) 


end 
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c c 

c AUTHOR : VIKAS LAKHOTIA c 

C DATE : FEB., 1990 c 

c C 

c THIS PROGRAM IS USED TO CALCULATE THE RESIDUAL STRESSES C 

C ARISING IN A CVD FIBER DUE TO THERMAL EXPANSION MISMATCH C 

C BETWEEN THE COATING AND THE CORE. THE TEMPERATURE USED IN C 

C THESE CALCULATIONS IS THE DIFFERENCE IN DEPOSITION C 

C TEMPERATURE AND ROOM TEMPERATURE. IT ACCOUNTS FOR THE C 

C ANISOTROPIC NATURE OF THE TiB2 COATING c 
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PROGRAM STRESS 

IMPLICIT DOUBLE PRECISION (A-H , O-Z) 

Open(ll ,File=' input . dat ’ , status= 'unknown* ) 

Ell, vll and etll are the Elastic modulus, Poissons ratio 
and Thermal expansion of tungsten. 

E12 , vl2 and etl2 are the Elastic modulus, Poissons ratio 
and Thermal expansion of titanium diboride parallel to c-axis. 
E22 , v22 and et22 are the Elastic modulus, Poissons ratio 
and Thermal expansion of titanium diboride parallel to a-axis. 
T is the deposition temperature . 

p is the ratio of fiber diameter to core diameter . 

Sa is the external stress applied, 
b is the final fiber diameter, 
r is the radial distance. 

Read(ll ,*) Ell, vll 
Read(ll,*) E12, vl2 
Read(ll ,*) E22, v22 
Read(ll,*) T, p, Sa 

et21=-0 . 205+5 . 959E-04+T+2 .77E-07*T**2-4. 369E-11*T**3 
et22=-0 . 118+3 . 092E-04+T+2 . 806E-07*T**2-3 . 996E-11*T**3 
etll=-4 . 3E-03+4 . 266E-04* (T-293)+8 ,479E-08*(T-293) **2 
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$ -1 .974E-ll*(T-293)**3 

C A negative sign is used in the folowing equations since 

C the fiber is being cooled. A factor of 0.01 is used as 

C the reported values of thermal expansion are in percantage. 

eti21=-0.010*(etll-et2l) 
etl22=-0.010*(etll-et22) 

AAl=Ell**2*E22*etl22 - Ell*E21*E22*etl22- Ell**2*E22*etl22*p**2+ 

- 2*Ell*E21*E22*etl22*p**2 - Ell*E21*E22*etl22*p**4 - 

- Ell*E21*E22*etl21*vl + 2*Ell*E21*E22*etl21*p**2*vl - 

- Ell*E21*E22*etl21*p**4*vl + Ell**2*E22*etl21*v21 - 

- Ell**2*E22*etl21*p**2*v21 + Sa* (- (Ell*E22*p**2*vl) + 

Ell*E22*p**4*vl + Ell*E22*p**2*v21 - Ell*E22*p**4*v2l) 
CCl=Ell**2*E21*etl21 - Ell*E21*E22*etl21 + 2*Ell*E21*E22*etl21* 

- p**2 - Ell**2*E21*etl21*p**4 - Ell*E21*E22*etl21*p**4 + 

- Ell*E21*E22*etl21*vl - 2*Ell*E21*E22*etl22*vl - 

- 2*Ell*E21*E22*etl21*p**2*vl + 4*Ell*E21*E22*etl22*p**2*vl + 

- Ell*E21*E22*etl21*p**4*vl - 2*Ell*E21*E22*etl22*p**4*vl + 

- 2*Ell**2*E21*etl22*v21 - 2*Ell**2*E21*etl22*p**2*v21 - 

- Ell**2*E2i*etl21*v22 + 2*Ell**2*E21*etl21*p**2*v22 - 

- Eil**2*E21*etl21*p**4*v22 + 

- Sa*(Ell**2*p**2 - Ell*E22*p**2 + Ell**2*p**4 + Ell*E22*p**4 + 

- Ell*E22*p**2*vl - EH*E22*p**4*vl + 2*Ell*E22*p**2*vl*v21 - 
2*Ell*E22*p**4*vl*v21 - 2*Ell**2*p**2*v21**2 - 
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- Ell**2*p**2*v22 + Ell**2*p**4*v22) 

DET=E11**2 - E11*E21 - E11*E22 + E21*E22 + Ell**2*p**2 + 

- Ell*E22*p**2 - 2*E21*E22*p**2 + Ell*E21*p**4 + E21*E22*p**4 + 

- Ell*E22*vl - E21*E22*vl - Ell*E22*p**2*vl + 

- 2*E21*E22*p**2*vl - E21*E22*p**4*vl - 

- 2*E21*E22*vl**2 + 4*E21*E22*p**2*vl**2 - 2*E21*E22*p**4*vl**2+ 

- 2*Ell*E21*vl*v21 + 2*Ell*E22*vl*v21 - 2*Ell*E21*p**2*vl*v21 - 

- 2*Ell*E22*p**2*vl*v21 - 2*Ell**2*v21**2 - Ell**2*v22 + 

- Ell*E21*v22 + Ell**2*p**2*v22 - 2*Ell*E21*p**2*v22 + 
Ell*E21*p**4*v22 

A1=AA1/DET 

C1=CC1/DET 

C2=C1/ (l-p**2)+Sa*p**2 
A2=A1/ (l-p**2) 

Sri=Al 

Stl=Al 

Szl=Cl 

Sz2=C2 

write (*,*) ’CORE STRESSES’ 
write(*,*) ’Radial Stress=Theta Stress’ 
write(*,10) Sri . 
write(*,*) ’Axial Stress=’ 
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write(*,*) Szl 

write (*,*) ’ COATING STRESSES’ 

write(* ,*) ’Radial Stress= (l-(b/r)**2)*’ 

write(*,10) A2 

write(*,*) ’Theta Stress= (l+(b/r)**2)*’ 
write(*,10) A2 

> 

write(*,*) ’Axial Stress=’ 
write(*,10) Sz2 
10 f ormat (2x , el6 . 8) 


end 
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Graphical representations 
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13.2 chromium Diboride Fibers 



In the deposition of chromium in the form of a thin film on an inert sub- 
strate many different precursors have been used.The main aim beeing to get 
the chromium in the gas phase at a relatively low temperature.The most 
important were the following. 

1. CrCl 2 1 

2. Cr 2 0 3 2 

3. Cr(CO) 6 3 , 


4. Cr 2 S 3 

Considering these and others the most appealing.considering availability and 
toxicity .were the following. 


Precursor | 

Melting point(C) | Boiling point (C) | 

CTO 2 CI 2 

-96.5 

117 

Cr{CO) 6 

83 

210 

: CrGh 

1150 

1300 


Chromyl Chloride as precursor 

In comparing the different compounds the use of chromyl chloride seemed 
very attractive because of its low melting point. As this compound had not 
been used in any previous research it was of major interest to us to see how 
it would behave combined with the carrier gas H 2 and the secondary reaction 
component BCl 3 .Hopefully it would be possible to deposit a pure film of CrB 2 
on our inert surface according to the overall equation : 

Cr0 2 Cl 2 + E 2 + BCl 3 = CtB 2 + ECl + - 


( 1 ) 




Samples with intermediate layer of carbon. 
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13*0. 0000000E+00, 

RHO 


1.000000 , 25*0.0000000E+00 

ERATIO 


F, 

OF 


F, 

FPCT 


F, 

FA 


F, 

MIX 

= 

20*0. 0000000E+00 , 

TP 

- 

T, 

HP 


F, 

s r 

= 

F, 

1 

d 

F, 
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F, 
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F, 
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so 

= 
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IONS 


F, 

I DEBUG 
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0, 

TRACE 
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SIUNIT 

vz 

F, 
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F. 

INHG 

SEND 

= 

F 
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NO INPT2 VALUE GIVEN FOR OF, EQRAT, FA, OR FPCT 
SPECIES BEING CONSIDERED IN THIS SYSTEM 
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J 3/65 
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J 6/61 

CLO 


J 3/61 

CL02 

J 9/65 

CL2 

J 12/65 

CL2 
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0 




T 
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J 6/73 

CR 

J 12/73 
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J12/73 
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CR03 
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H 
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HBO 
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J 
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J 3/64 

H02 

J 3/77 

H2 

H20 
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J 3/79 

H20(L) 





J 3/77 


J 

3/79 
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L 2/69 

H202 

J12/64 

H3B306 

0 




J 6/77 

OH 








J 

3/77 

02 


J 6/61 

03 

J 3/61 

H2(L) 



OF 

S 

0.000000 












EFFECTIVE 

FUEL 

EFFECTIVE OXIDANT 

HIXTURE 



ENTHALPY 

(KG-MOL) (DEG K) /KG 


HPP (2) 

-0.31158582E+03 


HPP (1) 

O.OOOOOOOOE+OO 


HSUBO 

-0.31158582E+03 


KG- ATOMS /KG 


BOP (I, 2) 

BOP (1,1) 

BO (I) 

CR 


0.10176423E-02 

O.OOOOOOOOE+OO 

0.10176423E-02 

CL 


0.1 72999 18E-01 

O.OOOOOOOOE+OO 

0.17299918E-01 

0 


0.20352845E-02 

O.OOOOOOOOE+OO 

0.20352845E-02 

B 


0.508&2113E-02 

O.OOOOOOOOE+OO 

0.50882113E-02 

H 


0.24423414E+00 

O.OOOOOOOOE+OO 

0. 24423414E+00 

1 20.735 

-29.635 

-53.604 -2.503 

-8.789 17.000 1000.000 


ADD CRB2 (S) 

1 -6.829 

-28.929 

-51.093 -5.720 

-8.802 8.000 1000.000 


ADD CRB (S) 

1 -11.116 

-29.355 

-52.810 -3.577 

-8.790 5.000 1000.000 


P *E CRB 2 (S) 

-10.546 

-29.222 

-52.374 -4.146 

-8.793 3.000 1000.000 


At>^B203 (L) 

V. -10.644 

-29.222 

-52.594 -4.048 

-8.793 3.000 1000.000 


2 -11.341 

-28.338 

-49.078 -2.937 

-8.915 4.000 1100.000 


ADD CRB2 (S) 

2 -10.740 

-28.199 

-48.677 -3.539 

-8.918 4.000 1100.000 


3 -10.452 

-27.298 

-45.171 -3.535 

-9.039 5.000 1200.000 


REMOVE CRB (S) 

3 -10.741 

-27.325 

-45.267 -3.391 

-9.038 3.000 1200.000 


4 -12.013 

-26.726 

-42.815 -2.666 

-9.147 5.000 1300.000 


5 11.604 

-39.326 

-91.239 -19.091 

-8.107 6.000 500.000 


ADD CRCL3 (S) 

5 -18.166 

-43.962 

-101.162 - • -4.206 

-8.068 6.000 '500.000 


REMOVE CRB2 (S) 

5 -27.888 

-40.721 

-94.569 -14.095 

-8.077 7.000 500.000 


ADD CRCL2 (S) 

5 -30.060 -39.997 

REMOVE CRCL3 (S) 

-93.113 -16.281 

-8.090 4.000 500.000 


5 -29.188 

-40.434 

-93.991 -14.963 

-8.081 4.000 500.000 




THERMODYNAMIC EQUILIBRIUM PROPERTIES AT ASSIGNED 
TEMPERATURE AND PRESSURE 


j vSE NO. 

1 



MOLES 

ENERGY STATE TEMP 

DENSITY 



CHEMICAL FORMULA 



;al/mol 
FUEL CR 

DEG K 
1.00000 CL 

G/CC 

2.00000 0 

2.00000 

1.000000 -1 

28656.680 
FUEL B 

G 300.00 

1.00000 CL 

0.0000 

3.00000 


5.000000 

96276.109 
FUEL H 

G 300.00 

2.00000 

0.0000 


120.000 

13.324 

G 300.00 

0.0000 



0/F“ 

0.0000 PERCENT FUEL- 

100.0000 EQUIVALENCE RATIO-12.2857 

PHI- 0.0000 

REACTANT 

DENSITY- 0. 

0000 




i 

THERMODYNAMIC PROPERTIES 


P, ATM 

1.0000 

1.0000 

1.0000 

1.0000 

1.0000 

T, DEG K 

1000.0 

1100.0 

1200.0 

1300.0 

500.0 

RHO, G/CC 

9.4172-5 

8.5462-5 

7.8167-5 

7.2142-5 1.9059-4 

H L/G 

-1.0036 

110.84 

227.11 

333.08 

-555.29 

1 ./G 

-258.16 

-172.53 

-82.706 

-2.6028 

-682.35 

G, v ^aL/ G 

-5547.52 

-6107.36 

-6677.89 

-7257.57 -2951.33 

Q . CAL/(G)(K) 

5.5465 

5.6529 

5.7542 

5.8390 

4.7921 

M, MOL WT 

7.727 

7.714 

7.697 

7.696 

7.820 

(dlv/dlp) t 

- 1.00000 

-1.00147 

-1.00002 

-1.00013 -1.00000 

(dlv/dlt)p 

1.0000 

1.0315 

1.0007 

1.0043 

1.0000 

CP. CAL/(G)(K) 

1.0056 

1.1941 

1 . 0408 

1.0861 

0.9721 

GAMMA (S) 

1.3436 

1.2955 

1.3305 

1.3152 

1.3540 

SON VEL, M/SEC 

1202.4 

1239.3 

1313.3 

1359.1 

848.4 


MOLE FRACTIONS 





CRB (S) 

7.7898-3 

4.3348-3 

0.0000 0 

0.0000 0 0.0000 0 

CRB2 (S) 

0.0000 0 

3.4403-3 

7.7535-3 

7.7544-3 0.0000 0 

CRCL2 (S) 

0.0000 0 

0.0000 0 

0.0000 0 

0.0000 0 7.8568-3 

BCL2 

4.219 -8 

3.553 -7 

1.749 -6 

9.684 -6 1.627-19 

BCL3 

1.7182-2 

1.3880-2 

1.0102-2 

9.8263-3 2.8373-2 

BH3 

1.337 -8 

7.277 -8 

2.284 -7 

1.123 -6 7.140-19 

BOCL 

9.030 -8 

1.026 -6 

7.338 -6 

4.168 -5 6.836-20 

B203 (L) 

1.6020-3 

1.7699-3 

2.3484-3 

2.0771-3 4.8026-3 

B303CL3 

3.5910-3 

3.4114-3 2.8053-3 

2.9959-3 .4,3532-4 

H 

2.129 -9 

2.416 -8 

1.834 -7 

1.028 -6 6.606-21 


3.993 -7 

4.888 -6 

‘3.568 -5 

2.344 -4 4.840-20 


7.654 -9 

9.899 -8 

8.866 -7 

4.546 -6 5.021-21 

'V 

7.0109-2 

8.0301-2 

9.3078-2 

9.3296-2 3.1427-2 

G 

Xj ’\ l 2o 

8.9973-1 

8.9286-1 

8.8387-1 

8.8375-1 9.2710-1 

8.727 -8 

3.301 -7 

1.164 -6 

2.234 -6 2.997-13 


ADDITIONAL PRODUCTS WHICH WERE CONSIDERED BUT WHOSE MOLE FRACTIONS WERE LESS THAN 0. 10000 
E-05 FOR ALL ASSIGNED CONDITIONS 
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Boron distribution among products 



H 2 /Cr0 2 Cl 2 













Chlorine distribution among products 













Chlorine distribution among products 
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Chromium iraction 


Chromium distribution among products 



BC1 3 / Cr0 2 Cl 2 
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Deposit v .ickness 


Deposit tTiiclcrxe s s vs. Time 



\ 1 ■ Linear dependence between deposition rate and Increase of the diameter of 

the deposited fiber 


T = 1 070°C 

H 2 /Cr0 2 Cl 2 = HO 
BC1 3 / Cr0 2 Cl 2 = 5.5 



) . 5 1.0 1-5 2.C 

Total gas flowrate (l/min) 


Deposition rate versus total gas flowrate 




BCl 3 /Cr0 2 Cl 2 


Stable solid phases at 1100°K 



Stable solid phases at 1100 K for the CVD system 





BV o /Cr0 2 Cl 2 


Kinetic experiments in CrB 2 region 



Stable solid phases at 1300 K with the points indicating at which conditions 
experiments were done. . 




BCl 3 /Cr0 2 Cl 2 








Proposed deposition rate from the second Model. 






Surface reaction verification 


O Untreated W-wire 
+ Treated W-wire 

1 hour <§> 1 00°C in HNO 3 (V0 -wt.%)- 


3.0 1.0 2.0 3.0 4.0 5.0 6.0 

Time (min) 


7.0 8.0 9.0 10.0 


The surface reaction was verified by using different surfaces 


Dubbel adsorp^on depositio 


Experimental vs. Dubbel adsorption 


model 


16.0 


15.0 


14.0 


I 13.0 



9.0 10.0 11.0 12.0 13.0 14.0 15.0 

Experimental deposition rate (m~ 6 /min) 


Agreement between the third model and the experimental results 



